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I discuss a Born (H 0) approximation of hadrons, motivated by a 
general feature of the data: The spectra of hadrons reflect their valence (qq 
or qqq) constituents, whereas hard scattering reveals a prominent sea quark 
distribution. Why do the sea quark d.o.f's not imply a richer spectrum? 
I look for an approach that can reconcile the quark and parton model 
descriptions of hadrons, and consider how this physics could emerge from 
the QCD Lagrangian. The possibilities are reduced by insisting that the 
approximation should be simple, yet adhere to the rules of quantum field 
theory. One might suspect that no such method exists - but the Born 
approximation presents itself. The description of relativistic bound states 
that it brings has interesting features which merit further exploration. 

PACS numbers: 12.38.Aw 



1. Introduction 

Hadrons have an intriguing double nature, which seems paradoxical 
based on our experience with non-relativistic bound states. On the one 
hand, meson and baryon quantum numbers are successfully classified in 
terms of their valence quark (qq, qqq) constituents [1]. Their gluon degrees 
of freedom are not evident, as shown by inconclusive searches for glueballs 
and hybrids. The quark model is generally successful in explaining masses 
and magnetic moments (with the notable exception of the pion mass). 

On the other hand, deep inelastic scattering shows that there are many 
gluons and sea quarks in hadrons. As seen in Fig. 1 there is an abundance 
of sea quarks even at the lowest Q 2 to which parton distributions may be 
reasonably extrapolated [2]. Since the proton is an ultra- relativistic bound 
state it is hardly surprising that it contains many light quark pairs. 

* Presented at the 50th Cracow School of Theoretical Physics, 9-19 June 2010 in Za- 
kopane, Poland. 
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Fig. 1. The Quark model (left) describes the hadron spectrum solely in terms 
of the valence constituents. The parton model (right) views hadrons as having 
important gluon and sea quark components. The sea quark and gluon distributions 
at Q 2 = 1 GeV 2 (lower right) are taken from [2]. 



So why do the sea quark degrees of freedom not contribute to the spec- 
trum? This is most likely related to their relativistic nature. Time-ordered 
perturbation theory shows that the intermediate states which contribute 
to scattering processes are frame-dependent. Fig. 2 illustrates how a co- 
variant Feynman diagram describing Coulomb scattering from an external 
source breaks into two time-ordered contributions. The relative weight of 
the second U Z" -diagram, which describes a particle pair fluctuation, de- 
pends on the frame. In the limiting case of the infinite momentum frame 
the Z-diagram is absent altogether. The Fock state structure of bound state 
wave functions must be similarly frame dependent. However, the quantum 
numbers of the bound state cannot change under boosts. Particle pair fluc- 
tuations apparently do not "count" as constituent degrees of freedom. 

The above argument is rather vague. It may, however, be formulated 
more precisely for electrons bound in a time-independent Coulomb potential 
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Fig. 2. Coulomb scattering from a time-independent source preserves the energy 
component p° = E of the scattering particle (covariant Feynman diagram on the 
left). In (t, p) space the amplitude is given by the sum of two time-ordered diagrams 
(right). The second "^''-diagram has a three-particle intermediate state. 
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Fig. 3. The Green function G formed by multiple interactions (K) of an electron 
propagating (S) in a static Coulomb potential. The initial and final electron mo- 
menta are denoted (p°, 0) and p = (p°,p). 

[3, 4]. The Green function G(p°,p) of the electron shown in Fig. 3 has poles 
at the bound state energies, say at p° = Er, 

G(p°,p) = S + SKS + SKSKS + . . . = S + SKG = + . . . (1) 

P u - E R 

The last equality displays this pole contribution, whose residue R(Er,p) 
satisfies the Dirac equation in momentum space. 

Due to the static nature of the potential the interactions transfer only 3- 
momentum, conserving the energy component p° of the electron momentum. 
Hence for p° > none of the intermediate electron propagators S(p°, k) 

depend on the ie prescription at p° = —E^, with = + m 2 . The 
Green function G(p°,p) in (1) is thus exactly the same whether a Feynman 
or Retarded prescription is used for the intermediate electron propagators, 

Sf ' r{p ' k) = \pO-E k + i s)( p o + E kTl e) (2) 



In particular, the energy Er of the bound state is unaffected by the ie 
prescription. 
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Expressing each intermediate propag ator S(p°,k) in (1) as a Fourier 
transform of the time-ordered S(t, k) reveals the intermediate states that 
contribute to the electron's bound state wave function at an instant of time 
t. Since Feynman propagators have support for t < they give rise to Z- 
diagrams as in Fig. 2. Retarded propagators on the other hand only allow 
forward motion in time, 

Sn(t, k) = M [(E fc7 ° - k ■ 7 + m)e-^ t + (E kl ° + k ■ 7 - rn)e iE ^] (3) 

and give no Z-diagrams. Hence there is no pair production and the bound 
state appears to contain a single electron, which can have either positive or 
negative energy. The distribution of this electron is given by the standard 
bound state wave functions determined by the Dirac equation. 

Even though this was a rather trivial example, the conclusion is inter- 
esting. The same bound state can be described using two quite different 
equal-time wave functions. The "true" picture is given by the Feynman 
propagator and leads to a wave function with Fock states containing any 
number of pairs. The momentum distributions of all those pairs should be 
specified to fully describe the bound state. On the other hand, the same 
bound state is described by the usual Dirac single particle wave function 
when retarded propagators are used. 

2. The h expansion 

The possibility to describe relativistic Dirac bound states using a single 
particle wave function, in spite of the multiple pairs they contain, is analo- 
gous the dual nature of hadrons depicted in Fig. 1. Specifically, the Dirac 
spectrum reflects only the degrees of freedom of the single electron. Can 
such a description be extended to bound states formed by the interactions 
of relativistic particles, without an external potential? 

The use of retarded propagators in the Dyson-Schwinger equation (1) 
was allowed due the absence of loops. All of the diagrams in Fig. 3 are 
tree-level, Born diagrams. Does the concept of "Born approximation" ex- 
ist for bound states? For scattering amplitudes the tree approximation is 
generally a good first approximation. Similarly, adding loop corrections to 
the propagators and vertices in Fig. 3 should not substantially change the 
nature of the Dirac bound states. 

Tree diagrams are generally regarded as the leading contribution to scat- 
tering amplitudes in the limit where the Planck constant h — > [5] . This is 
also the limit where one expects the laws of classical physics to apply. The 
fact that there are no classical bound states would suggest that there can 
be no Born level bound states, either. However, the relation between the 
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H — >■ limit and classical physics is not straightforward [6], as I shall next 
demonstrate for the harmonic oscillator [4]. 

The non-relativistic propagation of a particle from (ti,Xi) to (tf,Xf) in 



a potential V(x) 



2^,2 



i muj x 



is given by the path integral 



A(xi,xf,tf -U) = I [Dx(t)]exp 

[V£(t)] exp 



im 



dt{x 2 - u 2 x 2 ) 



2<t2\ 



(4) 



In the second equality the explicit dependence on ft was removed by scaling 
the coordinates as £ = x/y/h (the Jacobian is irrelevant for this discussion.). 
The full quantum mechanical bound state structure of the harmonic oscilla- 
tor thus persists as h — > in terms of the scaled variables £. The usual argu- 
ment that the rapid oscillation of exp(iS/h) selects classical paths for which 

the action S is stationary fails for propagation over distances of O , 

where the variation of the action itself is of O (h). It is only propagation 
between fixed (^-independent) positions Xi,Xf which becomes classical in 
the H — > limit, since that involves highly excited levels (E n = Huj(n + ^) 
with n of order 1/h). 

For the harmonic oscillator the concept of a Born approximation for 
bound states is rather trivial, since the dynamics in the rescaled variables is 
independent of h. In gauge field theory the h dependence is more interesting. 
Rescaling the fields of the QED action gives 



Sqed/H = I J d 4 x[^(i0-e4-m)^-lF^F^] 



d 4 x 



where in the latter expression 



^ = ip/Vh, A = A/^h. 



(5) 
(6) 



(7) 



The h dependence now appears exclusively in the interaction term, in the 
combination e\ffi. 

The dimensions of the fields and parameters in the action (5) can be 
easily worked out given that Sqed/^ is dimensionless, and the dimension 
of h in terms of energy E and length L is [h] = E ■ L (for c = 1). Thus 
[ij>] = E 1 / 2 ■ L-\ [A] = E 1 / 2 ■ L- 1 / 2 , [m] = 1/L and [e] = E^ 1 / 2 ■ L~ l l 2 . 

The coupling e appearing in the action S in (5) has a different dimension 
than the classical charge, [e] = E 1 / 2 ■ L 1 / 2 = [eh]. The dimensionless fine 
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structure constant is then 

e 2 e 2 h , . 

ol = - — - = - — 8 
4irh 47r 

Each loop correction brings a factor e 2 in the amplitude and thus also one 
power of h, if e is independent of H. Born (tree-level) approximations are 
thus obtained in the h — > limit at fixed e 1 . 




Fig. 4. Sum of ladder diagrams which reduces to the tree diagram sum of Fig. 3 
when the antifermion mass tends to infinity. 



The tree-level diagrams shown in Fig. 3, whose sum gives the Coulomb- 
Dirac bound states, stem in field theory from the crossed and uncrossed 
ladder diagrams shown in Fig. 4 [7]. One of the charged particles turns into 
the source of the Coulomb potential as its mass tends to infinity (in its rest 
frame). In contrast to what I said above, ladder diagrams with loops in this 
case give treedevel (Born) dynamics. This is possible since the bound state 
momenta depend on H. Atomic binding energies are of O (a 2 m) and hence 
ex H 2 according to (8). In order to stay at the position of a bound state pole 
as h — > the momenta must be reduced accordingly. Propagators with loop 
momenta of O {am) give inverse powers of momenta which cancel factors 
of a = e 2 h/4TT in the numerator. Hence all ladder diagrams contribute at 
the same order of a and h, allowing the perturbative sum to diverge at 
the bound state energies for any a. Loop corrections to propagators and 
vertices on the other hand give genuine higher order corrections in a and K. 



3. Hamiltonian formulation of Dirac states 

The QED Hamiltonian creates electron-positron pairs from the pertur- 
bative vacuum |0) in the presence of an external Coulomb potential. This 
reflects the Z-diagram contributions shown in Fig. 2, and implies that elec- 
tron bound states must have Fock components with an unlimited number 
of pairs. I shall now discuss how a Hamiltonian approach without pair pro- 
duction may be formulated, which corresponds to the retarded propagation 

1 In a study of how classical physics emerges from quantum field theory it may be 
more natural to keep the classical charge e fixed as h — > [6] . It then turns out that 
classical fields get contributions also from quantum loops. 
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discussed above. In the next Section I apply it to bound states formed by 
mutual interactions of relativistic particles [3] . Due to the absence of loop 
contributions this Hamiltonian formulation is equivalent to the standard 
one only at the Born level. 

Particle production is suppressed and retarded propagators obtained 
when the boundary condition is specified by the "retarded vacuum" , 

\0)R = N- 1 Hdl x \0) (9) 

p,A 

where the product is over all momenta p and helicities A of the antifermion 
creation operator. The normalization factor ./V is fixed by r(0|0}_r = 1. In 
the retarded vacuum 

bp,\\0)R = 4,aI°)« = and hence i/>(x)\0) R = (10) 

where ip{x) is the free (interaction picture) fermion field. Consequently the 
retarded propagator (3) is given by the standard operator matrix element 
in the retarded vacuum, 

S R {x-y)= R (Q\T[${x)${y)]\Q) R (11) 

The negative energy contribution to the propagator arises from the Sd 
term, which represents the removal of a positive energy antifermion from 
\0)r. The interaction Hamiltonian annihilates the retarded vacuum, 

Hi(t)\0) R = e J d 3 xA°(x)^(t,x)il>(t,x)\0) R = (12) 

which ensures the absence of particle production. 

A fermion bound state at t = may be parametrized in terms of its 
4-component Dirac (c-number) wave function <p(x) as 



\E,t = 0) = J d 3 x^(t = 0,x)<p(x)\0) R 

{2 ^3 P 2Ep E [ u] ^ X MP) b lx\°)R + vi (~P> X)p(p)d- P ,x\0) R ] (13) 

The negative energy components of tp(p) describe a state where cL p ,a has 
removed a positive energy antifermion from |0)_r. For the bound state to be 
stationary in time each Fock state amplitude (f>(t, x) must satisfy 

0(t, x) = R (0m, x)\E, t) = e- iEt ^(0, x) (14) 
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where 0(0, x) = <p(x) follows from f^ip a (t,x),ip^(t,x')^ = 5 3 (x — x') S a p 

and 4)(x)\0) R = 0. 

The time dependence of \E,t) is given by the interaction Hamiltonian 
(12). The stationarity requirement for the bound state at t = is then 

.d^M) = fl(0K ^M \ Ej o) + H ( |^(0, x)H T \E, 0) = EcP(0, x) (15) 

The interaction picture fields satisfy 

= 7 o ( _ zV ^ + mg) ^ (t) aj) (lg) 

and, making use of n{0\ip^(t, x) = 0, 

fl (0|V(0, 0) = e^VM*) (17) 

Using these relations in (15) gives the Dirac equation for the wave function 
ip(x) of a bound state of energy E in the external potential A°(x), 

(-iV • 7 + e7°AV) + m e )^(a;) = E-fip(x) (18) 



4. Hamiltonian formulation of mesons in QCD 

4-1. The importance of Coulomb exchange 

The above derivation of QED bound states in an external Coulomb po- 
tential, based on the retarded vacuum \0}r, gave the correct (Dirac) bound 
state energy and quantum numbers at the Born level. I now discuss how 
this method may be applied to Born level meson (qq) bound states in QCD. 
The wave function will describe a two-body state, whereas, like in the Dirac 
case, the same bound state formulated on the standard vacuum |0) has an 
indefinite number of sea quark pairs. 

The instantaneous Coulomb potential is determined by the positions 
of the charged constituents through Gauss' law. Since transverse gluons 
propagate in time their exchange involves intermediate states which include 
gluons. The derivation below concerns the valence qq Fock states only and 
thus applies insofar as Coulomb exchange is dominant. This is trivially 
the case in D = 1 + 1 dimensions, where there are no transverse degrees 
of freedom. However, the derivation applies also in D = 3 + 1 dimensions 
provided one allows for a homogeneous solution of Gauss' law (Section 4.4 
below) . This gives rise to a linear Coulomb potential which is of lower order 
in a s compared to the contribution of transverse gluon exchange. The linear 
potential is of course doubly welcome as a way of providing confinement. 
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4-2. The ud wave function x an d its gauge dependence 
The ud bound state at t = is analogously to (13) expressed as 

\E,t = 0} = J d 3 y 1 d 3 y 2 ^\t = 0,y 1 ) X AB (y 1 ,y 2 )^(t = 0,y 2 )\0} R (19) 

where now 

\0)r = N- 1 [J d^(p,\)bf(p,\)\Q) (20) 
p,X,A 

includes a product over colors. The wave function x which describes the 
state (19) is gauge dependent. Invariance of \E,t) under time-independent 
gauge transformations ip(t,x) — > U(x)ip(t,x) implies 

x(vi, y 2 ) -> u{y l )x{y l ,y 2 )u ] {y 2 ) (21) 

As an ansatz for a bound state solution I assume the existence of a gauge 
where the wave function is diagonal in color, 

x AB (yi,y2) = s AB x(yi,V2) (22) 

The stationarity condition for the bound state wave function corresponding 
to (14) is 

<t%jj>(t; Xl ,x 2 ) = R (0\i;$(t, x 2 )^g x (t, xi)\E, t) = e- iEt 4%${t = 0; Xl , x 2 ) 

(23) 

where </> a p(t = 0;xi,x 2 ) = Xap{xi,x 2 ). 

4-3. The qq matrix element of the gluon equations of motion 

The gluon field is determined by the equations of motion for each Fock 
state, i.e., for each qq configuration with the tt-quark of color C located at 
Xi and the ci-antiquark of color D at x 2 . The corresponding matrix element 
of the equation of motion, 

x + gfabcFPAl - g ^ ^fYT^f] \E,t) = (24) 

f=u,d 

should vanish to leading order in g. This gives 

X CD (x u x 2 )[d,Fr + gfabcFr A^] = g6 3 (x- x 1 )TC A 'y Yx AD (xi,X2) 

- g5 3 (x-x 2 )x CA (xi,X2h°J u T* D 

(25) 
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Only the color diagonal a = 3, 8 components of the Aa fields can be non- 
vanishing in the gauge where the wave function has the color structure (22). 
Then the commutator contributions vanish and the solution is abelian. 

In D = 1 + 1 the Coulomb potential ^3 8 is linear and we may set 
A\ 8 = 0. The bound state wave functions can then be derived as in Section 
4.5 below. 



4-4- The linear potential in 3+1 dimensions 

The field equations are consistent with an instantaneous linear potential 
even in D = 3 + 1 dimensions. The homogeneous solution 

A° a (x) = A 2 a £ a ■ x (o = 3,8) (26) 

trivially satisfies V 2 j4^ = for any (^-independent) values of the param- 
eters A a and unit vectors t a . Being of O (g ) this solution dominates the 
O (g) fields generated by the quark sources. Hence we may set A a = at 
leading order in g, ignoring transverse gluon exchange and the higher Fock 
components this would entail. Note that this is true also for relativistic 
constituents and in all coordinate frames. 

The fact that (26) (with A a = 0) solves the field equations (25) at O (g°) 
implies that the action is stationary under local variations of the fields. The 
action should also be stationary under variations of the global parameters 
A a , t a of the potential (26). This will ensure that the solution has rotational 
and color symmetry. 

The lagrangian density F^ v Fa V has an O (g°) contribution from the 
square of the homogeneous solution (26), an O (g) term from the interference 
of the homogeneous solution with the A® =3 8 fields generated by the quark 
sources in the field equations (25), and O (<? 2 ) contributions from squares 
of the O (g) fields. We need to consider the action at O (g) and thus also 
the field equations for ^4„ =3 8 including the O (g) quark sources, 

- V 2 A° a (x) = gT^ c [6 3 (x - x x ) - S 3 (x - x 2 )] (a = 3,8) (27) 

(no sum on the quark color C). The solution 

oT cc f 1 1 \ 
A° a (x; xi,x 2 , C) = A 2 a £ a -x + f - ] f (a = 3, 8) 

47T \\X — Xi\ \X — X2\ J 

(28) 

depends on the positions and color of the quarks, as well as on the param- 
eters A a , i a of the homogeneous solution. The instantaneous gluon action 
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with this Coulomb field is 



- 7 E / d'xF^Fr = \Y.l d * x ( VA °) 2 ( 29 ) 

a a 

= E I d " X + T a°ia ■ (*1 -X 2 )+0 tf) 

a=3,8 L J "'' 

The parameter Yla=3 8 A 4 , * s multiplied by the (infinite) volume of space. 
This term does not affect bound state evolution provided it is the same for 
all qq configurations. Hence we must require that 

A 4 = E A « ( 3 °) 

a=3,8 

is a universal constant, independent of X\,x 2 and the quark color C. The 
O (g) interference term in (29) is stationary wrt. variations of the unit vec- 
tors £ a provided i a || x\ — x 2 . Choosing i a = T^ c (xi - x 2 )/\T£°(xi - x 2 )\ 
gives (as seen below) an attractive linear potential. 

The (instantaneous) action (29) should be stationary also wrt. variations 
in the ratio A3/A8 which leaves A in (30) invariant. For quarks of color 
C = 1 the extremum of the O (g) term in (29) is obtained for A|/A| = y/S: 



Sin! = max j I (A3 + -j=Aij \x x - x 2 \j = ^j=\*l - x 2 \ (31) 

The calculation is the same for C = 2, whereas for C = 3 the extremum is 
obtained for Af/Af = 0: 



Sf n J z = max { ^\xi - x 2 \ \ = ^=\xi - x 2 \ (32) 



f pAf \ gA 2 

The factorization of the color dependence as in the ansatz (22) requires that 
the interactions are the same for all color components, i.e., that = 
Sgj 1 = S int . The equality of (31) and (32) was ensured by the color singlet 
nature of the action (29) and the color covariance of the equations of motion. 



4-5. Bound state equation 

Having determined the parameters A a and £ a in the A® potential (28) 
for each Fock state I proceed to impose the stationary time dependence (23) 
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on the bound state. Analogously to the Dirac case (15) I get 



<* CD (0;x 1 ,x 2 
-Jt = 



"a/3 VU ' a ' 1 ' JU2 ' _ viCD 



E^(0;x 1 ,x 2 ) = (33) 



+ R (0\^(0 1 x 2 )^ a (0,x 1 )[H I (0) - S mt ]\E,0) 

Here the Coulomb energy —Si n t stored in the gluon field contributes to the 
Hamiltonian with the opposite sign compared to the action (29). Using 
the time dependence (16) of the fields, the color structure (22) of the wave 
function and the interaction hamiltonian 

= 9 E / d 3 *i>f(t, x)A° a (x) T^ B i;f(t, x) (34) 

f=u,d 

gives the bound state equation 

7°(-iVi • 7 + m u )x(xi,x 2 ) - xi^i, x 2 )j°(iV 2 ■ 7 + m d ) 



= [E -V{x 1 ,x 2 )]x{xi,x 2 ) (35) 



Not surprisingly, it has the form of a "double Dirac equation" and as such 
was proposed a long time ago by Breit [8]. However, the approximations 
required to derive this equation from the underlying gauge theory have been 
clarified: lowest order in h (allowing the use of retarded propagation) as well 
as lowest order in the coupling g. The potential must be linear, 

V(xi,x 2 ) = ^^\xi-x 2 \ (36) 

where A is a free parameter with dimension of mass. In D = 1+1 dimensions 
even the perturbative potential is linear, with a coefficient oc g 2 . In D = 3+1 
the linear potential is obtained through a boundary condition for Gauss' law 
characterized by A through Eqs. (26) and (30). 



5. Discussion 

I have addressed an apparent paradox of hadron structure sketched in 
Fig. 1: The quantum numbers of hadrons reflect only their valence degrees 
of freedom (qq or qqq) even though direct DIS measurements indicate that 
there is a sea of quark pairs. This feature is closely related to the rela- 
tivistic nature of hadrons, and is seen also for an electron in an external 
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Coulomb potential. Dirac bound states have an infinite number of Fock 
components due to pair production in the perturbative vacuum. The same 
bound states are obtained at Born (no loop) level in a retarded vacuum, 
which suppresses pair production. The wave function then describes a sin- 
gle "valence" electron with both positive and negative energy components: 
The standard Dirac wave function. This example is encouraging as it shows 
that quantum numbers and mass spectra of relativistic bound states with 
an infinite number of constituents can be addressed analytically. 

Tree-level interactions are insensitive to the is prescription of propaga- 
tors, allowing the use of retarded boundary conditions. The basic approxi- 
mation may thus be formulated as being of lowest order in h, i.e., a Born 
approximation for bound states [3, 4]. This approximation is valid for rela- 
tivistic dynamics and provides a well defined starting point in the study of 
hadron structure. 

Another issue arises when the quark interactions are determined from 
the field equations rather than imposed by an external potential. Rela- 
tivistically moving charges generate transverse gauge fields which propagate 
with finite speed. Hence transverse gluon exchange involves intermediate 
Fock states with gluon constituents. This was demonstrated in QED for 
the Hydrogen atom [9] . Part of the binding energy of atoms in relativistic 
CM motion arises from transverse photon exchange and thus involves {epj) 
Fock states in addition to the \ep) state which dominates in the rest frame. 

A valence (qq or qqq) description of hadrons is thus only possible when 
the instantaneous Coulomb potential dominates transverse exchange. In 
D = 1 + 1 dimensions this is trivially true since there are no transverse gauge 
fields. In 3 + 1 dimensions the Coulomb field also dominates provided one 
allows for a homogeneous, O (<?°) solution of Gauss' law. Such a solution 
gives a constant color field extending to infinity and would seem to be 
unphysical. However, the instantaneous field is determined separately for 
each quark configuration and thus depends on the position and color of the 
constituents. The combined field strength far from color singlet hadrons 
in fact vanishes due to the coherent sum over quark colors. Translation 
invariance moreover requires that the hadrons be color singlets. 

The linear potential is of course very welcome also since it provides color 
confinement. I should emphasize that my analysis only shows that a linear 
solution is consistent with QCD at leading order. It does not provide a 
theoretical argument for why such a boundary condition must be used in 
QCD. To this day we also lack a proof that QED does not confine - which 
has not impeded theoretical progress. 

An analysis similar to the one outlined above for mesons can be carried 
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out also for baryons [3]. The potential turns out to be 

V(xi,x 2 , sea) = V (^i ~ x 2 ) 2 + ( x 2 ~ X3) 2 + (x 3 - xi) 2 (37) 

when the three quarks are located at x\,X2,x% and the color dependence 
of the wave function is given by e ABC . In the limit where two of the quarks 
are close together the potential (37) approaches the linear one (36) wrt. the 
position of the third quark. This construction is only possible when the 
number of colors equals the number of valence quarks. 

If the CM momentum k in the ud meson state (19) is separated as 

X (X 1 ,X 2 ) = e *(*l+*2)/2 Xk{xi _ X2) (38) 

the bound state equation reduces to 

- «V • [a, Xk(x)} + \k ■ {a, Xk{x)} + m u -f°Xk{x) - Xk{x)-f°m d 

= (E k - V)xk(x) (39) 

Because the bound state constituents are at equal time in all frames the 
fc-dependence of the wave function Xk is dynamical. However, the k- 
dependence of the bound state energy must be kinematic, Ej~ = Vk 2 + M 2 . 
Remarkably, the bound state equation (39) gives precisely this relation - and 
only for a linear potential [10]. The wave function Lorentz contracts in an 
intriguing way, since the boost parameter depends on the distance |x| be- 
tween the quarks, tanh £ = \k\/{E k — V{\x\)). Non-relativistic bound states 
{Ek 3> V) Lorentz contract in the standard way. 

The rest frame solutions of Eq. (39) have been studied using phenomeno- 
logical potentials [11, 12]. The radial wave functions can be singular at r = 
and at E — V{r) = 0. Solutions that are regular at both points have quan- 
tized masses which lie on asymptotically linear Regge trajectories. The wave 
function oscillates rapidly at large distances r between the quarks, where 
V{r) S> E. In this region the normalization per unit distance is constant, 
which may reflect the virtual pairs of the color string. 

The linear confining potential (36) appears as a "zeroth" order approx- 
imation in a perturbative expansion. Perturbative contributions of O {a s ) 
from gluon exchange, quark annihilation and loops still need to be added. 
This is consistent with the spirit of quark models, which typically assume a 
potential of the form 

V{r) = Cr-C F — (40) 
r 

where a non-perturbative linear potential is added to perturbative gluon 
exchange. This can be self-consistent provided a s is not too large even at 
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the confinement scale Q 2 = C, and if the linear potential may be effectively 
treated as a lowest order term in a perturbative expansion. 



Acknowledgements 

I am grateful and honored for the invitation to participate in this 50th 
Cracow School of Theoretical Physics. A long-time collaboration with Stan 
Brodsky on issues related to bound states is gratefully acknowledged, as 
well as helpful remarks by Stephane Peigne. 

REFERENCES 

[1] K. Nakamura et al. (Particle Data Group), J. Phys. G 37 (2010) 075021. 

[2] A. M. Cooper-Sarkar, Proc. 38th Int. Symp. on Multiparticle Dynamics 
ISMD08, Hamburg, Germany (2008), p. 449-453 [arXiv:0901.4001 [hep-ph]]. 

[3] P. Hoyer, arXiv:0909.3045 [hep-ph] and PoS EPS-HEP2009, 073 (2009). 

[4] S. J. Brodsky and P. Hoyer, arXiv:1009.2313 [hep-ph]. 

[5] J. Iliopoulos, C. Itzykson and A. Martin, Rev. Mod. Phys. 47 (1975) 165; 
C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill (1980), 
Sec. 6-2-1. 

[6] B. R. Holstein and J. F. Donoghuc, Phys. Rev. Lett. 93, 201602 (2004) [hep- 
th/0405239]. 

[7] D. R. Yennie, Brandeis Summer Institute Lectures on Strong and Electro- 
magnetic Interactions (Brandeis University, Waltham, Mass., 1963), Vol. I.; 
S. J. Brodsky, In *Brandeis Univ 1969, Proceedings, Atomic Physics and 
Astrophysics, Vol.1*, New York 1971, 91-169; and SLAC-PUB-1010 (1971); 
A. R. Neghabian and W. Glockle, Can. J. Phys. 61, 85 (1983). 

[8] G. Breit, Phys. Rev. 34 (1929) 553. 

[9] M. Jarvinen, Phys. Rev. D 71 (2005) 085006 [arXiv:hcp-ph/0411208]. 

[10] P. Hoyer, preprint HU-TFT-85-49 (1985) and Phys. Lett. B 172 (1986) 101. 

[11] D. A. Geffen and H. Suura, Phys. Rev. D 16 (1977) 3305. 

[12] W. Krolikowski, Acta Phys. Polon. B 23 (1992) 43. 



